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Abstract 

For a class of viable cosmological models in f{R) gravity in which deviation from the Einstein 
gravity decreases as a inverse power law of the Ricci scalar R for large R, an analytic solution for 
density perturbations in the matter component during the matter dominated stage is obtained in 
terms of hypergeometric functions. An analytical expression for the matter transfer function at 
scales much less than the present Hubble scale is also obtained. 
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I. INTRODUCTION 



Recent numerous observational data obtained from angular anisotropy and polarization of 
the cosmic microwave background radiation, large-scale gravitational clustering of galaxies 
and observations of distant supernovae explosions prove convincingly that the Universe ex- 
pands with acceleration at the present time while it was decelerating in the past for redshifts 
larger than about 0.7. If interpreted in terms of the Einstein general theory of relativity, 
this acceleration requires the existence of some new component in the right-hand side of 
the Einstein equations, dubbed dark energy (DE), which remains practically unclustered 
at all scales at which gravitational clustering of baryonic and dark non-baryonic matter is 
seen, and which effective pressure p^E is approximately equal to minus its effective energy 
density Ap_E- Thus, its properties are very close to those of a cosmological constant A (see 
Refs. l[ - 0] for some reviews). The simplest possibility of DE being exactly A com- 
bined with a non-relativistic non-baryonic dark matter (the standard spatially fiat ACDM 
cosmological model) provides a good fit to all existing observational dataQ]. In this case 
A acquires the status of a new fundamental physical constant. However, its required value 
is very small as compared to known atomic and elementary particle scales (not speaking 
about the Planck ones), so a firm theoretical prediction for this quantity from first princi- 
ples is lacking currently (although it may arise due to some non-perturbative effects, see e.g. 
Ref. [7|). 

On the other hand, in the second case when a component with qualitatively similar 
properties is assumed to exist - in the inflationary scenario of the early Universe, we are 
sure that this "primordial DE" may not be an exact cosmological constant since it should 
decay long time ago. That is why it is natural to assume by analogy that the present DE is 
not absolutely stable, too. 

An interesting alternative to the standard ACDM model is provided by "geometric" DE 
models based on f{R) gravity which modify and generalize the Einstein General Relativity 
by introducing a new function of Ricci scalar, f{R), into the gravitational field action instead 
of -R — 2A (see Ref. js] for a recent review). f{R) gravity, which in turn is a particular case 
of more general scalar-tensor gravity, contains an additional scalar degree of freedom or, in 
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quantum language, a massive scalar particle. ^ Because of this, viable models of present 
DE in f{R) gravity should satisfy a number of conditions which exclude many possible, in 
principle, forms of f{R). In particular, in order to have the correct Newtonian limit for 
R ^ R{to) ~ Hq where to is the present moment and Hq is the Hubble constant, as well 
as the standard matter-dominated FLRW stage with the scale factor behaviour a(t) oc t^/^ 
driven by cold dark matter and baryons, the following conditions should be fulfilled: 

\f{R)-R\^R, i?r(i?)<l, i?>i?o, (1) 

where the prime denotes the derivative with respect to the argument R. In addition, the 
stability condition f"{R) > has to be satisfied that guarantees that the standard matter- 
dominated FLRW stage remains an attractor with respect to an open set of neighboring 
isotropic cosmological solutions in f{R) gravity (in quantum language, this condition means 
that scalaron is not a tachyon).^ 

A number of functional forms has been proposed that can account for accelerated expan- 
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sion without A while passing laboratory and astronomical tests 

The present paper considers the evolution of density perturbation in f{R) model. In 
Section 2, we review the ACDM model. In Section 3, after discussing density perturbation 
in general f{R) gravity, we choose a specific inver se p ower-law form of f{R) — R which is 



the limiting form of the models proposed in Refs. 
solution for density perturbations. 
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ll| for R^ Rq and find an analytic 



II. DENSITY PERTURBATIONS IN THE ACDM MODEL 

Perturbed spatially fiat FLRW metric in the longitudinal gauge is 

ds"^ = -(1 + 2^)dt^ + a^it){l - 2^){dr'^ + r^dO + sin^ (2) 



^ This particle was dubbed "scalaron" in Ref. 9] where the particular case f{R) = R + /6M^ (plus small 
additional terms) was used to construct the first internally self-consistent inflationary scenario of the early 
Universe having a graceful exit to the subsequent radiation-dominated Friedmann-Lemaitre-Robertson- 
Walker (FLRW) stage through an intermediate matter-dominated reheating period. This inflationary 
model still remains viable since it predicts the slope of the primordial spectrum of scalar perturbations 
Hs and the tensor-to-scalar ratio r in agreement with the most recent observational data. 

^ The second stability condition f'{R) > which means that gravity is attractive and graviton is not a 
ghost is automatically fulfilled in this regime. 
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During the matter-dominated era, non-zero components of the energy-momentum tensor are 
given by 

T^ = -p-5p, Tf = -pd,v, (3) 

where p and 5p denote background matter density and its fluctuation, and v is the velocity 
potential of scalar perturbation. 

From the Einstein equations, we can derive the differential equation for comoving density 
perturbation defined by 

5 = ^ + mav. (4) 
Fourier transformation of density perturbation is given by 

where k denotes comoving wavenumber. In the following, we abbreviate 5k (t) just 6 for 
simplicity. The evolution of 6 in Fourier space is governed by 

6 + 2H6 - AnGpd = 0. (6) 

In matter dominated era, this equation becomes 

6 + -6 - ^6 = 0. (7) 

While it has two independent solutions proportional to t^/^ and t~^, we are only interested 
in the growing mode: 

2 

= 5ok (J-^ ' . (8) 

Here, ^ok is the initial value 5k (to)- 



III. DENSITY PERTURBATIONS IN f{R) GRAVITY 



We write the action in the following form: 



S 



167rG 



f{R) + C^"^^ 



(9) 



where G is the Newton constant and Z^^*^) is matter Lagrangian. If we take f{R) = R — 2A, 
Eq. (Q reduces to the Einstein-Hilbert action for the ACDM model. Below we consider 
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f{R) which vanishes for R = 0, so no cosmological constant is introduced by hand in flat 
space-time. 

In f{R) gravity, modified Einstein equations have the form: 



where 



F{R) 



dl_ 
dR' 



For dust-hke matter ([3]), the background equations take the form: 

?>FH^ = -{FR- f) -3HF + 87cGp, 
-2FH = F ~ HF + 8nGp, 
p + 3Hp = 0. 



(10) 

(11) 

(12) 
(13) 
(14) 



When deviation from the Einstein gravity is small, namely, f{R) = R and F{R) = 1, these 
equations yield the standard matter-dominated regime a{t) = ao(t/to)^^'^- 

The differential equation describing a density perturbation in the subhorizon regime is: 



5 + 2H5 - AnG^spS = 0, 



(15) 



where 



eff 



(16) 



The above equation was derived by Tsujikawa|l2|| (see also Ref. \L3i\) who took the subhorizon 
limit of the coupled equations of 6 and metric perturbations before reducing the system to 
a single equation.^ Its solutions were numerically studied in Refs. 15| - 18| in the linear 



regime and in Refs. 191] - 22|] in the non-linear one. 



On the other hand, de la Cruz-Dombriz et al. 23| derived the evolution equation of S 



without any approximation (see also Refs. 24, |25|). They then took the subhorizon limit to 
yield Eq. (33) in Ref. 23|. We find that, at the matter-dominated stage, their result can be 



In the two opposite limits ^ 1 and 1, Eq- (flSt reduces to the one derived in Ref. 14 1 

for matter density perturbations in general scalar-tensor cosmology for the values of the Brans-Dicke 
parameter ujbd ^ ^ and ujbd = correspondingly. 
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fc2 F,R ' ^ 



rewritten in the form of Eq. (fT5|) but with 

Apparently, Eq. ( fT6l) and Eq. ( ITTll have different dependence on fc. However, we can see 
that both of these two equations give practically the same result for observationally relevant 
f{R) models as follows. 

First, we consider two extreme cases with F^n/F ^ a'^/k'^ and F^r/F <^ a^/fc^. In the 
former case we see both equations have the same Gofr- In the latter case, Eq. ( IT6|) gives 
Geff = C/F, while Eq. ( TT71) reduces to Geff = G- However, since we are assuming the 
subhorizon regime, we find in this case F^r/F ^ a^/k"^ ^ H^^, so that F is practically 
constant (and in fact equal to unity) at the cosmological time scale. So, both equations give 
the same result in the two extreme regimes. 

Next we consider the intermediate regime with F^^/F ~ a'^/k'^. In this case, Eq. (fT6|) and 
Eq. f|T7|) do have different wavenumber dependence, but we still have F^r/F ~ /k"^ ^ 
and again F = 1 in cosmological time scale. It has been numerically confirmed in Ref. [23i] 
that in such observationally viable f{R) models density fluctuations evolve in the same way 
whichever equations one uses. Hence both equations give the same result in this case, too. 

We therefore continue to use Eq. (fT5|) with Eq. (fT6l) . From now on, we adopt a specific 
f{R) model such that 

F{R)^f\R) = l-[f^\ , n>-l (18) 



with i?o ~ Hi (but still _Ro < Riio))- This f{R) corresponds to the models in [lO|, lUl in 



the regime R^ Rq, and we shall use it in this regime only.^ Equation (flS!) then becomes 

U±S-^^'-±^^^I^J = 0. (19) 

3t 3t2 1 + 3A(t/to)2'^+^/3 



with 



A{n,k)=^^^^±^. (20) 



The parameter n used here has the same sense as in Ref. [lo| . while it is equal to 2n in the notation of 
Ref. [ll|. 
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A. Asymptotic behavior 



5k(t) = 5ok(-) . (21) 



5k(t) = So^C{k) {^] . (23) 



We can read asymptotic behavior from the differential equation Eq. (fTOll . 
(i) t^O 

In this hmit, we can neglect y4(t/to)^"^*''^ respect to 1. Therefore, Eq. f|T9|) reduces to 
Eq. ([7]). This reproduces the result of the ACDM model. The solution is Eq. (jH]), 

2 

t \ 3 



(ii) t ^ oo 

Then Eq. f|T9|) reduces 

In this regime the growing mode is given by 

-1+733 

The coefficient C{k) is the transfer function for matter perturbations which will be 
derived from the analytic solution found below. 



B. Analytic solution 

By changing the variable from t to r = (t/to)° with a = 2n + 8/3, Eq. (HM can be 
rewritten as 

Here a prime denotes derivative respect to r. This equation can be solved in terms of the 
hypergeometric function 2-Fi- The two independent solutions are 

- (a - 1: , (1) ") . 

In the following discussion, we consider the upper sign case only, because the other solution 
corresponds to the decaying mode of perturbations and is singular at t — )■ 0. 
Let us check the asymptotic behavior of the solution, Eq. ( 125|) . 
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(i) t^O 



^ 5ok ( ^ 

^0 



(26) 



(ii) t — )■ oo 



5k(t) 5ok X 



r(i + |^)r(^ 



r 1 + 



5+V33 
6a 



5+V33 
6a 



3(n + 1)A;2 
an-Ro 



-5+^33 
6cK 



-1 + ^33 
6 



(27) 



In both limits, the asymptotic behavior agrees with that one given by Eq. ([21]) and 
Eq. ( 123|) . respectively. Furthermore, here we can read off the transfer function, C{k), which 
appears in Eq. fl23l) : 



Cik) 



r (i + 2(3It4)) r ( 



2(3n+4) 



r 1 + 



5+V33 
4(3n+4) 



5+V33 
4(3n+4) 



3(n+ l)fc^ 



-5+\/33 
4(3n+4) 



(28) 



IV. CONCLUSIONS AND DISCUSSION 



We have obtained an analytic solution describing the growth of density perturbation at 
the matter-dominated stage for a specific class of viable cosmological models in f{R) gravity. 
Initially, the solution behaves in the same way as in the ACDM model, while it experiences 
an anomalous growth at late times (redshifts of the order of a few). We also find an analytic 
expression for the matter transfer function which shows that an initial perturbation power 
spectrum acquires the additional power-law factor oc A;^"* with 

-5 + ^33 



An., 



3ra + 4 



(29) 



at scales much less than the present Hubble scale, as originally shown in Ref. Hl].^ Clearly, 
this additional factor is absent in the matter power spectrum at the recombination time. So, 
by comparing the form of the primordial matter power spectrum derived from CMB data 
and from galaxy surveys separately, it is possible to obtain an important constraint on the 
parameter n characterizing this class of cosmological models in f{R) gravity, although we 
do not have much stringent constraints on it at present 261]. If we take an upper limit on 



^ Note the difference in the notation for n between Ref. 



11| and our paper which was mentioned above. 
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Aus as An^'^^ = 0.05, which is a conservative bound for now ll|, and assume that Rq is not 
much less than (if otherwise Rq <^ H^, deviation of the background FLRW model from 
the ACDM one is very small), we obtain a constraint 



n > 4.96 



Anl 



0.05 



- 1.33. 



(30) 



Future observational data together with a more detailed theoretical analysis may well yield 
a more stringent bound on n. 

Of course, the f{R) gravity model f|T8l) is viable for a finite range of R only, in particular, 
for R^ Rq. For R ~ -R(to) ~ H^, it has to be substituted by a more complicated expression 
admitting a stable (or, at least metastable) de Sitter solution, e.g. by the models presented 



in Refs. 
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III]. As a result, the equation for matter density perturbations has to be solved 



numerically for recent redshifts z < 1. However, evolution in this region may add only a 
/c-independent factor to the total matter transfer function. Therefore, the k exponent in 
Eqs. (!28|) and (!29|) does not depend on a concrete form of f{R) for R ~ -R(to)- 

Also, the model f|T8|) should not be used for too large values of R for several reasons. 
First, the effective scalaron mass squared = 1/3Fr (in the regime F^rH^ <^ 1) grows 
quickly with R to the past and may even exceed the Planck mass making copious production 
of primordial black holes possible. As was noted in Ref. this problem may be avoided 
by adding the i?^/6M^ term to f{R) that bounds the scalaron mass from above just by 
M. Simultaneously, such a change of f{R) at large R removes the "Big Boost" singularity 
(in terminology of Ref. 27| where such a singularity appeared in a different context), which 
generic appearance in the model ( fTSi) was shown in Ref. 28|]. However, the value of M 
should be sufficiently large in order not to destroy the standard cosmology of the present 
and early Universe. In particular, its values considered in Refs. 29|, l30l] seem not to be high 
enough for this purpose. Indeed, the simplest way to solve one more problem of this f{R) 
cosmological scenario (also noted in Ref. |ll]) - overproduction of scalarons in the early 
Universe - is, as usual, to have an inflationary stage preceding the radiation-dominated one. 
Then M should not be smaller than H during its last 60 e-foldsjand for M ^ 3 x 10^'^ GeV 
the scalaron will be the infiaton itself according to the scenario 91] . 

Note that for sufficiently large n the corresponding correction to F{R) may become more 
important than the second term in the right-hand side of Eq. (fT8|) already at the matter- 
dominated stage. For example, even for M as large as 3 x 10^^ GeV, the term R/3M'^ 



becomes larger than {Rq/R)"'~^^ for R = 3x 10^''-R(to) (corresponding to the matter-radiation 
equahty) if n > 10. However, this does not affect the exact solution obtained in the previous 
section since at this moment F^^/ F ^ o? /k"^ for all scales of interest, so Ges = G in Eq. (IT^ 
irrespective of an actual structure of F{R) — 1. In turn, if F^^/F > a'^/k'^, the R^ correction 
is negligible for all scales of interest at the matter-dominated stage. Therefore, this high- 
R correction to the model (ITSjl needed to obtain a viable cosmological model of the early 
Universe does not change our results. 
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